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ABSTRACT 

A quasi-three-dimensional thin-layer Navier- 
Stokes analysis was used to predict heat transfer 
to rough surfaces. Comparisons are made between 
predicted and experimental heat transfer for tur- 
bine blades and flat plates of known roughness. 
The effect of surface roughness on heat transfer 
was modeled using a mixing length approach. The 
effect of near-wall grid spacing and convergence 
criteria on the accuracy of the heat transfer pre- 
dictions are examined. An eddy viscosity mixing 
length model having an inner and outer layer was 
used. A discussion of the appropriate model for 
the crossover between the inner and outer layers 
is included. The analytic results are compared 
with experimental data for both flat plates and tur- 
bine blade geometries. Comparisons between pre- 
dicted and experimental heat transfer showed that 
a modeling roughness effects using a modified mix- 
ing length approach results in good predictions of 
the trends in heat transfer due to roughness. 

N omencla t ure 


Pr t - Turbulent Prandlt number 
p - Pitch between roughness elements 
Re - Unit Reynolds number 

5 - Reference surface area 
Sf - Total frontal area 

s - Surface distance 

St - Stanton number based on inlet conditions 
U - Inviscid velocity 
u r - Friction velocity 
y - Distance normal to surface 

- Normalized distance, {yU t /v)\/C / /2 

6 - Full boundary layer thickness 

k - Von Karman’s constant 

A - Roughness density parameter 

p - Dynamic viscosity 

p t - Turbulent eddy viscosity 

v - Kinematic viscosity 

v t - Turbulent eddy viscosity, (Mt/p) 

p - Density 

oj - Vorticity 


Af - Frontal area 
A B - Windward wetted area 
A+- Damping coefficient 
Cj - Friction factor 

D - Diameter of hemispherical roughness 
h - Roughness height 
/i eq - Equivalent roughness height 
L - Distance between roughness elements 


t - Length scale 


Subscripts 

c - Edge of boundary layer 
m - Measured 
p - Predicted 


This paper is declared a work of the U.S. Government and is 
not subject to copyright protection in the United States. 
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INTRODUCTION 

The effect of inaccuracies in predicting turbine 
blade temperatures significantly affects turbine 
durability. As discussed by Graham* 1 ), a 50°C 
underprediction of blade temperature can decrease 
service life by a factor of two. Taylor* 2 ) measured 
the surface roughness of in-service turbine blades. 
His results showed that for typical Reynolds num- 
bers the surfaces are not hydraulically smooth. 
Since actual turbine blades are rough, heat transfer 
predictions assuming smooth surfaces can signifi- 
cantly underpredict external heat transfer. This 
in turn could result in underprediction of the 
blade surface temperatures. Both two-dimensional 
boundary layer analyses, such as the STANS code 
of Crawford and Kays* 3 ), and Navier-Stokes anal- 
yses, such as those of Chima* 4 ) and Hah* 5 ), have 
been used to predict smooth turbine blade heat 
transfer. Models to account for surface roughness 
in boundary layer analyses have been developed 
by Cebeci and Chang* 6 ), and by Taylor, Coleman, 
and Hodge* 7 ). Navier-Stokes analyses have the ad- 
vantage over boundary layer analyses in that they 
can be used to predict blade heat transfer beyond 
the point of separation. The present work dis- 
cusses the results obtained when a surface rough- 
ness model was incorporated into the quasi-three- 
dimensional thin-layer Navier-Stokes analysis de- 
veloped by Chima* 4 ). 

There are two different philosophies with re- 
gard to modeling surface roughness. The approach 
taken by Taylor et al. * 7 ) explicitly accounts for 
roughness geometry. This approach accounts for 
blockage due to the roughness elements, the in- 
creased drag due to their presence in the flow field, 
as well as a correlation to account for additional 
heat transfer between the fluid and the discrete 
roughness elements. Consequently, this approach 
requires detailed knowledge of the surface rough- 
ness, as well as empirical correlations for the drag 
and heat transfer to the roughness elements. Hosni, 
Coleman, and Taylor* 8 ) have used this model in a 
boundary layer analysis to predict heat transfer for 
a rough flat plate over a wide range of Reynolds 
numbers. The other philosophy, as expressed by 
Cebeci and Chang* 6 ), is to model the roughness 
effect based on the offset in the near-wall velocity 
profile. In this model the mixing length is modified 


to account for the effective height of the roughness. 
The effective height is a function of both roughness 
density and shape. Computationally, the model of 
Cebeci and Chang is significantly easier to imple- 
ment into a Navier-Stokes analysis than the more 
detailed model of Taylor et al. Also, Hosni et al.* 8 ) 
reported that their boundary layer calculations, us- 
ing the model of Taylor et al., were done using 250 
points within the boundary layer. This is far more 
than is practical to use in a Navier-Stokes analysis 
for turbine blade heat transfer. Consequently, the 
model of Cebeci and Chang was implemented in 
the Navier-Stokes analysis. 

Navier-Stokes heat transfer predictions can be 
affected by the near-wall grid spacing. The sen- 
sitivity of the heat transfer predictions for both 
smooth and rough surfaces will be examined. Be- 
cause roughness promotes very early transition, 
heat transfer to rough surfaces can be assumed 
to be fully turbulent. The effect of different ap- 
proaches to turbulent eddy viscosity modeling in 
the leading edge stagnation region is discussed. 
Comparisons are made with experimental data to 
demonstrate the suitability of the roughness mod- 
els for predicting turbine blade heat transfer. The 
data of Hosni, Coleman, and Taylor* 8 ) for flow over 
a roughened flat plate was used for comparisons. 
Comparisons were made for the heat transfer on 
turbine blades using both the data of Liebert* 9 ) 
and the data of Blair and Anderson* 10 ). 

METHOD of ANALYSIS 

The roughness model of Cebeci and Chang* 6 ) 
was incorporated into the Navier-Stokes analysis of 
Chima* 4 ) by modifying a version of the Baldwin- 
Lomax * n ) turbulent eddy viscosity model. In this 
model there is an inner and outer region. The tur- 
bulent eddy viscosity in the inner region is given 
by: 

(Mt)iNNER = fit* M (!) 

The length scale, £, is given by: 

t = K y{ l-e"* + ' A+ ) (2) 

Where y + = yu/u T and A + = 26. In the outer 
region the turbulent eddy viscosity is given by: 

(/*t)oUTER = C'outerP^’wake^kleb (3) 
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Where Cquter = 0.0269, and Fwake is given by: 


fwAKE = ™*«(ltaAX^MAX>-25yMAXVDIF/*MAx) 

( 4 ) 

The term «dif is the maximum velocity along 
the near normal grid line out from the surface. The 
quantities j/max an d i'MAX are found from: 

F(y) = »H(1 - (5) 

Vmax is the value of y for which F(y) is a maxi- 
mum. Fkleb is given by: 

^kleb = [l + 0.004(y/yMAx) 6 ] (0) 

For rough surfaces Cebeci and Chang use an 
offset distance, Ay, which is added to the physi- 
cal normal distance. This was implemented by re- 
placing y by y + Ay in the above equations. For 
example, equation 2 for a rough surface becomes: 

1=k { y + Ay) (l - e -(v + +*v + )M + ) (7) 

At the wall y = 0, but t > 0 when the surface is 
rough. Therefore, > 0 at the wall. The heat 
transfer predictions were made using an effective 
wall thermal conductivity. Pr t was assumed equal 
to 0.9 for the determination of the effective waill 
conductivity. 

Cebeci and Chang give Ay as: 

Ay + = 0.9 (yfht q - k+ e‘ fc V«) (8) 

Cebeci and Chang state that equation 8 is valid in 
the range 4.535 < h+ q < 2000. Figure 1 shows the 
variation of Ay+ with h+ q . When the equivalent 
roughness height, h equals 4.535, Ay + is zero. 
This automatically accounts for the fact that small 
roughnesses, < 5), are hydraulically smooth. 
For large values of h^ q the slope of the curve pro- 
gressively decreases. Consequently, the effect of in- 
creasing roughness height becomes less significant 
as the roughness becomes large. 

In the boundary layer analysis of smooth walls 
the length scale, £, is held constant when ny ex- 
ceeds 0.0865. In the Navier-Stokes analysis the 
boundary layer thickness, 5, is not clearly defined, 
and is not used in the Baldwin-Lomax turbulence 


model. For rough surfaces the turbulent eddy vis- 
cosity in the outer region was calculated in two 
ways. First, y was unmodified in the outer-region 
eddy viscosity equations. In the second way, y 
was replaced by y + Ay in the outer-region eddy 
viscosity equations. The consequence of these as- 
sumptions on the rough surface heat transfer will 
be discussed. 

The Baldwin-Lomax turbulence model has an 
abrupt crossover model. In this model, when 
(Mt)oUTER > (pt) INNER, Mt is taken as the outer 
value. Granville^ 12 ) discusses two other crossover 
models which blend the inner and outer regions. 
These two models are: 


Vj_ 

V 



, (*v /D inner 
tann . . r 

OUTER (^r/^JoUTER 


(9a) 
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The effect on the predicted heat transfer using each 
of these models will be discussed. 

RESULTS and DISCUSSION 
Heat transfer predictions are known to be sen- 
sitive to the near-wall grid spacing in a Navier- 
Stokes analysis. This sensitivity could be differ- 
ent depending on whether the surfaces are rough 
or smooth. Figure 2 shows the sensitivity of heat 
transfer predictions to the near-wall spacing for a 
representative rotor geometry for both smooth and 
rough surfaces. The results are given in terms of 
a reference y~*". The same definition is used as was 
used by Boyle in reference 13. The definition is: 


!/ref = 0.17yi/2e°' 9 / s 0 ' 1 (10) 


Where yi is the distance from the surface of the 
first grid line. Re is the exit Reynolds number per 
unit of length. y£ EF is only a weak function of 
the surface distance, $, and s is taken as the axial 
chord. The same definition is used for the reference 
surface roughness, EF except that yi is replaced 
by h tq . 


*ref = o.m eq Re°7s° 1 (11) 

The calculations are for fully turbulent flow. 
The smooth wall results shown in figure 2a are not 
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very sensitive to the near wall spacing. A greater 
sensitivity was shown in reference 13 for a similar 
comparison. The sensitivity was greater when the 
flow was laminar, or close to transition. 

Figure 2 shows a somewhat greater sensitiv- 
ity to the near wall spacing for the rough surface 
than for the smooth one. The reference values were 
calculated assuming a smooth surface, and the cal- 
culation depends on the inverse of the friction fac- 
tor. Thus, for rough surfaces the reference values 
is lower than one based on the correct friction fac- 
tor by the ratio >/(C/) rough f{C /) smoo th * The 
rough surface friction is not known before the anal- 
ysis is done, and equation 11 is useful for determin- 
ing the grid spacing. Consequently, the rough wall 
heat transfer shown in figure 2b is more sensitive 
to the near-wall spacing than the smooth wall data 
shown in figure 2a. 

The primary purpose of this investigation was 
to examine the effects of surface roughness on heat 
transfer predictions. Consequently, grids with close 
near-wall spacing were used. Calculations were 
typically done using 145 x 54 C-grids. Approxi- 
mately 5000 iterations were needed to insure con- 
vergence. The computations took about 700 CPU 
seconds on a Cray-XMP. It is expected that further 
optimization of the code and calculation procedure 
would result in faster convergence. 

Figure 3 shows the effect of the different 
crossover models for fully turbulent flow. The cal- 
culations are for a smooth wall. The difference in 
the calculated heat transfer among the different 
models is not great. Of the two blending mod- 
els, equation 9b is closer to the Baldwin- Lomax 
crossover model. The Baldwin- Lomax model re- 
sults in a heat transfer prediction that has greater 
fluctuations between adjacent nodes than either of 
the other crossover models. This is consistent with 
the abrupt nature of the crossover model. For some 
cases the solution stability or convergence might be 
enhanced using a blended crossover model. 

The equivalent roughness height, h tqy is a 
function of the actual roughness and the rough- 
ness density. A number of correlations have been 
proposed to obtain h cq . Sigal and Danberg^ 14 ) dis- 
cuss various correlations for obtaining the equiv- 
alent height. The ratio of equivalent to actual 
roughness height, h eq /h y is correlated as a func- 
tion of the roughness density parameter, A. Fig- 


ure 4 shows their results with different correlations 
for both two- and three-dimensioned roughness. A 
is defined differently for the different correlations. 
The definition used by Sigal and Danberg was: 

Asb = ^(t:) " i 12 “> 

The definition used by Dvorak ( 15 ) as interpreted 
by Simpson^ 16 ) was: 

Adv = jr ( 126 ) 

The definition used by Dirling^ 17 ^ was: 

=!(£)"” < 12 ‘) 

Sigal and Danberg reasoned that three-dimensional 

roughness has a lower h eq because flow is able to go 
around as well as over a three-dimensional obstacle. 
Since different definitions were used for A, Dvorak’s 
correlation does not necessarily result in a lower 
value of h cq . 

In the correlation of Sigal and Danberg, and in 
the correlation of Dvorak, for the same roughness 
shape, A is proportional to the height-to-pitch ratio 
for two-dimensional roughness. However, for three- 
dimensional roughness A is proportional to the 
square of this ratio. Consequently, even for closely 
spaced roughness, only two-dimensional roughness 
results in A values in the region where the slope of 
the correlation is positive. 

Hosni et al.^ 8 ^ conducted experiments on 
smooth and rough flat plates. The rough surface 
test sections were formed from machined aluminum 
plate, and had integral hemispherical roughness el- 
ements. They were of diameter D, and were spaced 
one of three distances, L y apart to form a three- 
dimensional roughness array. This resulted in three 
different roughness densities. Table I shows the 
equivalent height ratio for different spacings using 
the three definitions of A. Because A was defined 
differently by the different investigators, there is 
less variation in the equivalent height ratio than 
might be expected from figure 4. The equivalent 
height ratio varied from about 0.02 to somewhat 
greater than 1.0 as LjD varied from 2 to 10. 
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Figure 5 shows comparisons of predicted and 
measured heat transfer over a flat plate for differ- 
ent equivalent roughness heights. Since the Navier- 
Stokes analysis used here was developed for turbo- 
machinery applications, the flat plate results were 
calculated as a thin uncambered airfoil in a low 
solidity cascade. The airfoil had an elliptical lead- 
ing edge with a 10:1 radius ratio. This was done 
to avoid any separation in the leading edge region, 
as might occur with a circular leading edge. The 
experimental data are shown for clarity as indi- 
vidual points, but were derived from a curve fit 
through data. As the Reynolds number increases 
the predicted heat transfer becomes greater than 
the experimental data for the smooth surface. The 
predicted heat transfer is greater by nearly 10%. 
In the Navier-Stokes analysis the solidity may not 
have been low enough. There was some small ac- 
celeration of the freestream flow due to boundary 
layer growth. This would increase the heat transfer 
above that for an isolated flat plate. 

Table II compares the ratio of predicted to 
measured heat transfer for the different roughness 
spacings at a Reynolds number of 400000. Com- 
parisons are given for the values of h eq /h from 
the different roughness correlations. Also shown 
in this table is the ratio of predicted to measured 
heat transfer for the smooth surface. The analy- 
sis predicts the trends in the experimental data. 
Generally, the best agreement is found when the 
lowest value of h cq /h for a given spacing is used 
in the analysis. The most recent correlation of 
roughness data is that of Sigal and Dan berg. This 
is also the correlation that agrees best with the 
data. Sigal and Danberg showed that spherical 
segment data might deviate significantly from their 
proposed correlation. This is the roughness geom- 
etry that most closely resembles the roughness el- 
ements used by Hosni et. al. A line using the 
extremes of the spherical segment data would be 
about a factor of 3 less in equivalent height ratio 
than the three-dimensional correlation of Sigal and 
Danberg. Consequently, it is not suprising that a 
heat transfer analysis using the correlation of Sigal 
and Danberg would overpredict the heat transfer 
when hemispherical roughness elements are used. 

There is little rough surface turbine blade 
heat transfer data available that is suitable for 
comparisons with analytic predictions. Suitable 


data would have well defined information on the 
roughness characteristics, and not just the rough- 
ness height. Measurement techniques suitable for 
smooth blades may not be suitable for rough sur- 
faces. For example, Taylor, Taylor, Hosni, and 
Coleman^ 18 ) showed that surface mounted sensors 
which do not have the same roughness character- 
istics as the blade itself may indicate heat transfer 
rates different from that of the blade. If the surface 
of the blade is artificially roughened using a low 
conductivity material such as sand, this material 
may effectively move transition forward on the sur- 
face, but also form a thermal resistance. These re- 
quirements suggest that experimental heat transfer 
data for rough turbine blades should include base- 
line smooth surface data. The experimental data 
of Tarada* 19) showed that for many cases there 
was no increase in heat transfer for rough blades 
when the flow is turbulent. In these data, in which 
roughness was achieved by coating a smooth blade 
with sand, the roughness only tripped the bound- 
ary layer close to the leading edge. The early trip- 
ping of the boundary layer resulted in an increased 
average surface heat transfer. However, for many 
cases the heat transfer was lower in the fully turbu- 
lent region for rough blades than for smooth blades. 
The analysis does not predict this behavior. The 
predicted effect of roughness for turbine blades is 
similar to the predictions for a flat plate shown in 
figure 5. 

In contrast to the data of Tarada, the data 
of Blair and Anderson* 10) show higher heat trans- 
fer for a roughened blade in the fully turbulent re- 
gion. While the roughness height was given, no 
information was given to determine the roughness 
density. Sand was used to roughen the surfaces, 
but the tests were conducted in a much different 
thermal environment. The significance of an insu- 
lating layer may have been less in this environment 
than in the tests of Tarada. The tested blade was 
the rotor of a large-scale low-speed annular cas- 
cade. The turbulence intensity is high in front of 
the blade due to the presence of the upstream sta- 
tor. Figure 6 compares predicted and measured 
heat transfer for the smooth surface. Two predic- 
tions are shown in figure 6a. The one in best agree- 
ment with the data was for a turbulence intensity 
of 10%. The augmentation of laminar heat trans- 
fer was included until the flow was fully turbulent. 
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The transition model, which is very sensitive to the 
turbulence intensity, is the one described in refer- 
ence 13. The other calculation neglects the effect 
of turbulence with respect to both transition and 
augmentation of laminar heat transfer. Compar- 
isons between the two calculations show the effect 
of freestream turbulence to be large. 

Rough surface predictions are shown in figure 
6b. Calculations are shown for different values of 
the equivalent height ratio. Since sand was used to 
roughen the blade, an equivalent height ratio near 
unity is expected. In fact this does give reasonably 
good agreement with the data. The pressure sur- 
face heat transfer calculation is less sensitive than 
the suction surface to the assumed value of equiv- 
alent height ratio. 

Figure 7 shows heat transfer predictions for 
both smooth and rough surfaces of the turbine ro- 
tor blade tested by Liebert^. Only a rough sur- 
face blade was tested. The primary purpose of the 
test was to verify a heat flux measurement tech- 
nique. Measurements were made on a 1:1 scale 
SSME high-pressure fuel turbine rotor blade. Be- 
cause of the relatively small size only one sensor 
was installed on each blade. Only three blades were 
instrumented in the test. To obtain a complete 
mapping of the blade heat transfer, several blades 
could be instrumented. The surface roughness pro- 
files were measured. Both the surface roughness 
height, A, and the roughness density, A, were cal- 
culated. The equivalent height ratio, h eq /h , was 
determined to be between 0.2 and 0.4. Predictions 
are shown for different height ratios as well as for a 
smooth blade. The tests were conducted in a tur- 
bine tester rig. This rig matched the engine turbine 
inlet temperature, but operated at a lower pres- 
sure than the actual engine. The lower pressure 
resulted in lower than engine Reynolds number. If 
the roughness height were the same, the value of 
in the engine would be higher them in these 
tests. These tests, like those of Blair and Ander- 
son, were also in a high turbulence environment. 
Because of the high turbulence intensity and high 
Reynolds number, fully turbulent predictions are 
shown. The data are in reasonably good agree- 
ment with the prediction for an equivalent height 
ratio of 0.1. The figure also shows that had the 
roughness been of the same height, but of a dif- 
ferent density, the heat transfer might have been 


significantly higher. 

The fully turbulent predictions shown in fig- 
ure 7 used an averaging process near the stagnation 
point. There is symmetry in flow properties on ei- 
ther side of the grid line that intercepts the stagna- 
tion point. Consequently, the calculated vorticity 
is nearly zero along this line. This in turn results 
in a turbulent eddy viscosity, /**, that is nearly zero 
along the grid line which leads up to the stagna- 
tion point. For grid lines adjacent to this line, the 
vorticity is large, as is /**. In the analysis /it was 
taken as the average of the values at five grid lines 
in the stagnation region. 

CONCLUSIONS 

Predictions for rough surface heat transfer ex- 
hibit a greater sensitivity to the near wall spacing 
than do those for smooth surfaces. This results in 
the need for closer near wall grid spacing for the 
rough surface predictions. 

The modified mixing length approach of Ce- 
beci and Chang can be used to predict the trends in 
heat transfer rates with roughness. The disagree- 
ment between the analysis and the experimental 
data could generally be accounted for by the un- 
certainty in the assumptions regarding roughness. 
The primary uncertainty is the equivalent height. 

The results of this investigation illustrate the 
importance of having a well characterized surface 
definition for experiments used to validate rough 
turbine blade heat transfer predictions. Not only 
should the roughness height and density be known, 
but the thermal characteristics of the roughness 
should also be known. For the data to be suit- 
able for comparison with analytic predictions it is 
necessary to quantify any additional thermal resis- 
tivity caused by the roughness. 
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Table I. - Equivalent height ratios for roughness geometries of Hosni et al. 



Equivalent Height Correlation 


Sigal &Danberg 

Dvorak 

Dirling 

LjD 

A 

w* 

A 


m 

h eq /h 


30.8 

1.58 

10.2 

2.43 

10.0 

|iHS9 

m 

123 

0.257 

40.7 

0.50 

KUB 

Blfa 

H 

772 

0.023 

257 

0.061 

■ 

EH 

0.079 


Table II. - Ratio of predicted to measured heat transfer for a 
flat plate at Re x = 400000. Data of Hosni et al. 



Equivalent Height Correlation 


Sigal &Danberg 

Dvorak 

Dirling 

LjD 

^eq/^ 

st p /st m 

^eq/^ 

st p /st m 

h e q/h 

st p /st m 

2 

1.58 

1.64 

2.43 

1.85 

1.69 

1.67 

4 


1.36 

0.50 

1.57 

0.45 

1.52 

10 

0.023 

1.01 

0.061 

1.15 

0.079 

1.22 

Smooth 


1.10 


1.10 


1.10 



Fig. 1 Offset in mixing length as a function 
of roughness height. 
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STANTON NO. X 1000 



-1 0 1 
SURFACE DISTANCE, s/s T 

a) Smooth surface 




Fig. 3 Effect of crossover model on blade heat transfer 
for smooth blade 
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b) Rough surface, h* REF =60. 

Fig. 2 Effect of near-wall spacing on calculated 
blade heat transfer distribution. 


Fig. 4 Equivalent roughness height ratio as a 
function of roughness density parameter 
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Fig. 5 Calculated and measured heat transfer for 
a smooth and rough flat plate 
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Fig. 6 Comparison of predicted and measured rotor blade 
heat transfer. Data of Blair and Anderson. 
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Fig. 7 Comparison of predicted and measured rotor blade 
heat transfer. Data of Liebert. 
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